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Abstract. High resolution spectra of 34 halo population dwarf and subgiant stars have been obtained with 
VLT/UVES and used to derive sulphur abundances from the AA8694.0, 8694.6 and AA9212. 9, 9237.5 Si lines. In 
addition, iron abundances have been determined from 19 Fell lines and zinc abundances from the AA4722.2, 4810.5 
lines. The abundances are based on a classical ID, LTE model atmosphere analysis, but effects of 3D hydrody- 
namical modelling on the [S/Fe], [Zn/Fe] and [S/Zn] ratios are shown to be small. We find that most halo stars 
with metallicities in the range —3.2 < [Fe/H] < —0.8 have a near-constant [S/Fe] ~ +0.3; a least square fit to 
[S/Fe] vs. [Fe/H] shows a slope of only —0.04 ± 0.01. Among halo stars with -1.2 < [Fe/H] < —0.8 the majority 
have [S/Fe] ~ +0.3, but two stars (previously shown to have low a/Fe ratios) have [S/Fe] ~ 0.0. For disk stars 
with [Fe/H] > — 1, [S/Fe] decreases with increasing [Fe/H]. Hence, sulphur behaves like other typical a-capture 
elements, Mg, Si and Ca. Zinc, on the other hand, traces iron over three orders of magnitude in [Fe/H], although 
there is some evidence for a small systematic Zn overabundance ([Zn/Fe] ~ +0.1) among metal-poor disk stars 
and for halo stars with [Fe/H] < —2.0. Recent measurements of S and Zn in ten damped Lya systems (DLAs) 
with redshifts between 1.9 and 3.4 and zinc abundances in the range —2.1 < [Zn/H] < —0.15 show an offset 
relative to the [S/Zn] - [Zn/H] relation in Galactic stars. Possible reasons for this offset are discussed, including 
low and intermittent star formation rates in DLAs. 

Key words. Stars: abundances - Stars: atmospheres - Galaxy: evolution - Galaxies: high-redshift - quasars: 
absorption lines 



1. Introduction 

Sulphur is generally regarded as an a-capture element. 
The work on Galactic stars by Frangois (1987, 1988) sup- 
ported this view by showing that [S/Fe] increases from 
zero at solar metallicities to a plateau level of about +0.5 
dex in the metaUicity range —1.8 < [Fe/H] < —1.0. This 
is an analogous behaviour to those of other a-elements, 
Mg, Si, and Ca - see Norris et al. H2001|l and Carretta et 
al. H2002|l . who both find [a/Fe] to be nearly constant at 
a level of +0.4 dex in the metaUicity range —4 < [Fe/H] < 
— I and then to decrease towards the solar [a/Fe] ratio 
for [Fe/H] > — 1. The standard interpretation is that this 
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trend arises from the time delay in the production of about 
two thirds of the iron by supernovae (SNe) of Type la rel- 
ative to the near-instantaneous release of the a-elements 
by Type H SNe. In this connection it should be noted 
that the metaUicity at which Type la SNe start to con- 
tribute (the "knee" of the [a/Fe] - [Fe/H] trend) seems to 
depend on the orbital properties of the stars. Halo stars 
with [Fe/H] ~ —1 belonging to the outer halo tend to 
have lower [a/Fe] than stars in the inner halo (Nissen & 
Schuster [TW7I Stephens & Boesgaard HUU^ Gratton et 
al. 12003(1 . Stephens & Boesgaard find a slope of —0.15 for 
the mean [a/Fe] vs. [Fe/H] for 53 stars in the metaUic- 
ity range —3.8 < [Fc/H] < —0.6, but this is for a sample 
with special kinematics, i.e. stars with large maximum dis- 
tances from the Galactic center and/or the Galactic plane 
or stars having extreme retrograde motion. 
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Recent observations of the AA8694.0, 8694.6 A Si lines 
in spectra of metal-poor stars by Israelian & Rebolo ( 200l| 
have, however, challenged the view that sulphur is an 
a-element. Their data, obtained with the 4-m William 
Herschel Telescope on La Palma, suggest that [S/Fe] in- 
creases linearly with decreasing [Fe/H] to a level as high 
as [S/Fe] ~ +0.7 at [Fe/H] = -2.0. The study of Takada- 
Hidai et al. 1)2002(1 , based on Keck High Resolution Echelle 
Spectrograph (HIRES) observations, supports a quasi lin- 
ear dependence of [S/Fe] on [Fe/H] although in their case 
[S/Fe] reaches only +0.5 dex at [Fe/H] = -2.0. 

As a possible explanation of the high value of [S/Fe] 
in metal-poor stars, Israelian & Rebolo (2001 ) proposed 
that very massive supernovae with exploding He-cores and 
a high explosion energy make a significant contribution 
to the early chemical evolution of galaxies. According to 
Nakamura et al. ()2001jl these hypernovae overproduce S 
with respect to O, Mg and Fe. With this intriguing pos- 
sibility in mind a more thorough investigation of sulphur 
abundances in halo stars seems worthwhile. 

A clarification of the trend of S abundances is also 
much needed in deciphering the chemical enrichment of 
damped Lya systems (DLAs), widely regarded as the 
progenitors of present-day galaxies at high redshift. The 
importance of sulphur stems from the fact that, un- 
like most other heavy elements, S is not depleted onto 
dust. Consequently, observations of the relatively weak S ii 
triplet resonance lines at AA1250, 1253, 1259A yield a di- 
rect measurement of the abundance of S in DLAs. Another 
element for which this is the case is Zn, and indeed most 
of our current knowledge of the chemical evolution of the 
universe at high redshift is based on surveys of [Zn/H] in 
DLAs (e.g. Pettini et aL lTMni Prochaska & Wolfe 2002). 
// S is an a-capture element, then its abundance relative 
to Zn (assumed to be an iron-peak element) could be used 
as 'a chemical clock' to date the star-formation process at 
high z. Specifically, if a major star formation episode in 
a DLA occurred within « 0.5 Gyr prior to the time when 
we observe the galaxy, we would expect to measure an 
enhanced [S/Zn] ratio, and vice versa. Measurements of 
[S/Zn] in DLAs have been relatively scarce until recently, 
but are now becoming available in increasing numbers 
thanks the Ultraviolet and Visual Echelle Spectrograph 
(UVES) on the ESQ Very Large Telescope (VLT). The 
VLT/UVES combination affords high spectral resolution 
and high efficiency over most of the optical spectrum, from 
blue and near-UV wavelengths to the far red. However, 
without a secure knowledge of the behaviour of S in metal- 
poor Galactic stars we clearly stand little chance of inter- 
preting the situation at high z. 

In the present paper we report on a survey of sul- 
phur abundances in 34 metal-poor dwarf stars based on 
high resolution observations of the AA8694.0, 8694.6 and 
AA9212.9,9237.5A Si lines obtained with UVES. The 
AA4722.2,4810.5A Zni lines are also included in our spec- 
tra, allowing us to study the Galactic evolution of both 
sulphur and zinc. 
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Fig. 1. Portion of the VLT/UVES spectrum of HD 110621 
{V = 9.9, [Fe/H] = -1.66) in the far red. Top panel: 
the spectrum before removal of the numerous telluric H2O 
lines which occur in this region. Bottom panel: the spec- 
trum after division by that of the B-type star HR5488 
using the IRAF task telluric. The Si triplet, a Mgii 
line and the Paschen C Hi line of HD 110621 can all be 
clearly seen. Note that the spectrum shown has not been 
corrected for the radial velocity of the star (which intro- 
duces a shift of -t-6.4 A), and that the wings of the Hi line 
have been removed by the continuum fitting. Excessive 
noise in the stellar spectrum is seen where strong H2O 
lines have been removed. 

2. Observations and data reduction 

2.1. UVES observations 

Previous determinations of sulphur abundances in metal- 
poor stars, including those by Israelian & Rebolo (,,2001,1 
and Takada-Hidai et al. H2002|l . have been based on the 
high excitation (xoxc = 7.87 eV) AA8694.0, 8694.6 S i lines. 
However, these lines are very weak in metal-poor dwarfs 
and giants (with equivalent widths W < S mA at [Fe/H] ~ 
—2.0), and become vanishingly small at [Fe/H] ^ —2.5. In 
our survey we have taken advantage of the high efficiency 
in the far red of the VLT/UVES instrument (Dekker et al. 
I2000f) and concentrated on the stronger S i triplet (xcxc = 
6.52 eV) at AA9212.9, 9228.1, 9237.5 A. The dichroic mode 
of UVES was used to cover the spectral region 6700 — 
10000 A in the red arm and 3750 — 5000 A in the blue arm; 
the latter region includes a number of Fe 11 lines suitable 
for determining the iron abundance, as well as the Zni 
lines at 4722.2 and 4810.5 A. 

The observations were carried out using the UVES im- 
age sheer #1 (Dekker et al. I2002|l which has an entrance 
aperture of 2.1 x 2.7 arcsec and makes 3 slices along the 0.7 
arcsec wide entrance slit with a total length of 8.0 arcsec. 
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Fig. 2. Left: A sequence of spectra around the AA8694.0, 8694.6 Si lines. Right: A sequence of spectra around the 
A9212.9 Si line. The arrows indicate wavelengths where telluric H2O lines have been removed. 



The resulting resolution of the spectra is A/ A A ~ 60 000 
with 4 pixels per spectral resolution element. 

Program stars were selected from the Stromgren pho- 
tometric catalogue of Schuster & Nissen p988p supple- 
mented with a few very metal-poor stars from Ryan et al. 
lUnSSIl- The selection criteria were: 5700 < T^s < 6500 K, 
3.7 < log (7 < 4.4 and a smooth distribution of metallici- 
ties from [Fe/H] = -3.2 to -0.6. The UVES spectra were 
obtained in service mode during the period March - June, 
2001. The total integration time for a V = 11 mag star 
was about 60 min, split into three separate exposures so 
that cosmic ray hits could be removed by comparison of 
the three spectra. The sky background could be checked 
adjacent to the stellar spectrum, but was insignificant (ex- 
cept for emission lines) even for the faintest star (G64-12, 
V = 11.46). Typical S/N ratios are 300 in the blue spec- 
tral region, 250-300 at the AA8694.0, 8694.6 S i Uncs, and 
150-200 in the region of the AA9212.9, 9228.1, 9237.5 Si 
triplet. 

2.2. Reduction of the spectra 

The blue spectra were reduced with standard MIDAS rou- 
tines for order definition, background subtraction, flat- 
field correction, order extraction and wavelength calibra- 



tion. Bias, dark current and scattered light corrections are 
included in the background subtraction. The spectra were 
then normalized by a continuum function determined by 
fitting a spline curve to a set of pre-selected continuum 
windows estimated from the solar atlas. Finally, correc- 
tion for the radial velocity of the star, measured from 20- 
30 narrow lines, was applied before the measurement of 
equivalent widths. For a few stars, it was checked that re- 
duction with NO AO's IRAF echelle package gave nearly 
identical results. 

The red spectra were reduced in the same way as the 
blue spectra but using the IRAF echelle package. A prob- 
lem with the 9212-9238 A region is the presence of nu- 
merous, strong telluric H2O lines, see Fig. ^ To remove 
these lines, fast rotating, early-type stars were observed 
on each night and reduced in the same way as the pro- 
gram stars. The IRAF task telluric was then applied 
to remove the telluric lines. When the early-type star is 
observed at about the same airmass as the program star 
this technique also serves to remove a residual fringing of 
a few percent in the red region (A > 7000 A) remaining 
after the flatfielding of the spectra. 

A sequence of spectra of stars spanning the whole 
metallicity range of the sample is shown in Fig. El As 
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seen, the Si line at 9212.9 A can be clearly detected in a 
[Fe/H] 3.2 star. 



2.3. Equivalent widths 

Equivalent widths of the Fe II, Zn I and S i lines were mea- 
sured by Gaussian fitting, or direct integration if the fit 
was poor. For the present set of data, we could estimate 
empirically the accuracy of our equivalent widths by com- 
paring two independent sets of measurements for nine 
stars (HD 103723, HD 104004, HD 110621, HD 121004, 
HD 140283, HD146296, G 16-13, G 66-30, and G64-12) 
which were observed twice (on different nights). The first 
set of observations was carried out on March 9 and 10 
while the second set was made from March 12 to May 9, 
2001. The total exposure time for each star was the same 
at each epoch. The comparison of the Fe ii and Zn i lines 
is shown in Fig. |3| As can be seen from the figure, the 
data fall close to the 1:1 line with a standard deviation of 
only ±0.84 mA. This corresponds to an error of ±0.6 mA 
for one measurement of an equivalent width in the blue 
part of the UVES spectra. 

The weak Si lines at 8694.0 and 8694.6 A could be 
measured in the more metal-rich part of our sample 
([Fe/H] > -1.5). Among the stronger triplet lines A9228.1 
falls very close to the center of the Paschen C Hi line 
and its equivalent width could therefore not be measured 
in a reliable way. The other two S i lines are, however, 
ideal for abundance determination in the metallicity range 
-3 < [Fe/H] < -1. 

The comparison of equivalent widths of S i lines 
measured from spectra obtained on different nights is 
shown in Fig. 0] Again, the data are evenly distributed 
around the 1:1 line, but with a larger scatter for the 
AA9212.9, 9237.5 lines (cr = ±3.3 mA) than is the case for 
the AA8694.0, 8694.6 pair (cr = ±0.8mA). The reason is 
the lower S/N in the 9212 ~ 9238 A region, made worse 
by the telluric lines. Thus, the equivalent widths of the 
AA9212.9, 9237.5 lines have been measured with a preci- 
sion of ±2.3 niA (one spectrum), whereas the precision for 
the weaker AA8694.0, 8694 pair is as good as ±0.6 mA. 

The error of the measurement of the equivalent widths 
of the AA9212.9, 9237.5 S i lines is particular important for 
estimating the precision of the sulphur abundance deter- 
minations for the most metal-poor stars, where the lines 
are weak. An independent check of the estimated error 
may be obtained by comparing the sulphur abundances 
derived from the A9212.9 and the A9237.5 Si line, respec- 
tively. For nine stars with [Fe/H] < — 2 having equivalent 
widths of the S i lines ranging from 5 to 36 niA the mean 
difference of the sulphur abundances derived from the two 
sets of lines is 0.03 dex with a rms scatter of the differ- 
ence of ±0.08 dex, only. This is, in fact, a smaller disper- 
sion than expected from the estimated equivalent error of 
±2.3 mA, which should then be considered as a conserva- 
tive estimate. 
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Fig. 3. A comparison of equivalent widths of Fell and 
Zn I lines measured from two sets of spectra for nine stars 
observed on March 9-10 and March 12 - May 9, 2001, 
respectively. 

The measured equivalent widths for 19 Fe ii lines, the 
two Zn I lines and the four S i lines in 34 program stars 
are listed in Table Al^. When no value is given it is either 
because the line is too strong or too weak to provide a 
reliable abundance, or, in the case of the AA9212. 9, 9237.5 
S I lines, is affected by residuals from the removal of strong 
telluric H2O lines. In addition to the stars listed in Table 
Al, three stars, HD 99382, BD -13 3834 and G 18-54, were 
observed, but they turned out to be double-lined spec- 
troscopic binaries, and are excluded from our abundance 
analysis. 

3. Stellar parameters 

The determinations of stellar parameters follows the 
method described by Nissen et al. I|2()()2|l . Here we briefly 
outline the procedure and discuss the precision of Teg, 
log (7 and [Fe/H]. The final values of Tea and logg are given 
in Table 121 together with values for [Fe/H] and the micro- 
turbulence as derived from Fell lines fSect. ITT|l . As the 
cahbrations of Tea and logg depend somewhat on [Fe/H], 
the determination of the parameters is an iterative pro- 
cess. 



Table Al is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
,http: / /cdsweb.u-strasbg.fr /cgi-bin/qcat? J / A-|- A /..■/, .. 
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Fig. 4. A comparison of equivalent widths of the 
AA8694.0, 8694.6 and AA9212.9, 9237.5 A Si lines mea- 
sured from two sets of spectra for nine stars. 



3.1. Effective temperature 

Toff was determined from the b — y and V ~ K colour in- 
dices using the IRFM calibrations of Alonso et al. H1996|l 
as modified by Nissen et al. H2002|l . The Stromgren uvby-(3 
photometry was taken from Schuster & Nissen (,1988,) for 
the large majority of the stars supplemented with unpub- 
lished photometry of Schuster et al. (2003) for the remain- 
ing stars. The (6 — y)o calibration of Schuster & Nissen 
H1989|l including a zero-point correction of +0.005 mag 
(Nissen'1994') was used to derive the interstellar reddening 
excess, li E{b — y) > 0.015, the reddening is considered sig- 
nificant and the V magnitude as well as the mi and the 
ci indices are corrected according to the relations given 
by Nissen et al. (2002). These values and the reddening 
estimates are given in Table 

The K photometry was taken from Carney ()1983|) . 
Alonso et al. (|1994fl and the Two Micron All Sky Survey 
(2MASS)^. In this connection we note that the 2MASS K 



^ This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by 
NASA and the National Science Foundation. K magnitudes 
were taken from the Second Incremental Release. Additional 
values from the recent full release have only a marginal effect 
on the derived T^s and abundances. 



magnitudes are on the so-called iiT-short system (Cutri et 
al. lM7S|l . whereas those of Carney (1983) are on the CIT 
(California Institute of Technology) system and Alonso et 
al.'s (,1994,) magnitudes are on the TCS (Telescope Carlos 
Sanchez) system. Small differences between these systems 
may exist. Carpenter (2001) has derived an average trans- 
formation X(2MASS) = K(CIT) - 0.024 without any 
colour term. For our actual sample of metal-poor turnoff 
and subgiant stars the agreement is even better. Twelve 
stars with both CIT and 2MASS photometry have a mean 
difference in K of 0.002 only with a rms-dispersion of the 
difference of 0.026. For 16 stars with both Alonso et al. 
and 2MASS photometry the mean K difference (TCS - 
2MASS) is -0.019 with a dispersion of 0.030. Hence, the 
three systems agree within ±0.02 mag in K for our sam- 
ple, which according to the Alonso et al. (|1996fl Toff cal- 
ibration oi V — K corresponds to an error of ±20 K in 
Tcff. The corresponding effect on the derived abundances 
is small compared to other error sources as seen from 
Table 01 We have, therefore, adopted the straight mean 
of the K magnitudes if two or three sources were avail- 
able. Table ^ lists the corresponding V — K value cor- 
rected for interstellar reddening according to the relation 
E{V-K) = 2.7 E{B-V) = 3.8E{b-y) (Savage & Mathis 
IMS), if E{b-y) > 0.015. 

Comparing the Tcs values determined from b ~ y and 
I^ — ii' we find a mean difference (A = Tcs{b—y)—Tcs{V — 
K)) < A >= +44 K for 25 stars with a standard deviation 
of (t(A) = 78K (one star). Similar values were found for 
the sample of metal-poor stars in Nissen et al. f2 002fl . 
The mean value of Tcs{b — y) and Tcs{V — K) is adopted 
for the 25 stars. For nine stars without K photometry, a 
value Teff(& — j/) — 22 K has been adopted. The typical 
observational errors are 0.007 mag in 6 — y and 0.05 mag 
vnV — K^ which correspond to an error of ±50 K in Teg 
in either case. Taking into account the uncertainty in the 
reddening estimate (also corresponding to about ±50 K 
in Teff), we estimate the 1-a statistical error of Tcff to be 
around 70 K. 

3.2. Surface gravity 

The surface gravity was determined from the fundamental 
relation 

M Teff , 



log- 



50 



log 
0.4(Mfc„i 



4 log 



T 



ff,0 



(1) 



bol,Q) 



where M is the mass of the star and Mboi the absolute 
bolometric magnitude. 

The absolute visual magnitude My was determined 
from a new calibration of the Stromgren indices derived by 
Schuster et al. (j20()3|l on the basis of Hipparcos parallaxes, 
and also directly from the Hipparcos parallax (ESA [TM7|l 
if available with an error (j{'k)/'k < 0.3. Column 12 and 13 
of Tablc^list the photometric and the parallax based val- 
ues of My. They compare reasonably well for the 18 stars 
having My par. Excluding one star, G 186-26, with a large 



6 



P.E. Nissen et al.: Sulphur and zinc abundances in Galactic stars and damped Lya systems 



Table 1. Stromgren photometry, colour excess, V — K index, T^s from b — y and V — K, spectroscopic value of [Fe/H] , 
and absolute magnitudes derived from the Stromgren photometry and from the Hipparcos parallax including an error 
corresponding to the parallax error. If E{b — y) > 0.015, the V magnitudes and the photometric indices have been 
corrected for interstellar absorption. 
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.582 


-0.005 


1 


.448 


5822 


5824 


-0.73 


4.20 


4.22 


0.08 


HD 160617 


8.733 


0.347 


0.051 


0.331 


2 


.584 


0.011 


1 


.408 


5953 


5909 


-1.79 


3.27 


3.42 


0.31 


HD 179626 


9.210 


0.373 


0.095 


0.293 


2 


.588 


0.015 


1 


.560 


5750 


5649 


-1.14 


4.24 


3.59 


0.39 


HD 181743 


9.687 


0.351 


0.052 


0.224 


2 


.582 


0.000 






5885 




-1.93 


4.97 


4.95 


0.34 


HD 188031 


10.148 


0.328 


0.058 


0.310 


2 


.592 


-0.001 






6076 




-1.79 


4.02 






HD 193901 


8.660 


0.383 


0.099 


0.221 


2 


.568 


-0.002 


1 


.533 


5652 


5691 


-1.12 


5.16 


5.46 


0.12 


HD 194598 


8.354 


0.344 


0.091 


0.269 


2 


.588 


-0.011 






5928 




-1.17 


4.63 


4.62 


0.15 


HD 215801 


10.044 


0.334 


0.049 


0.330 


2 


.572 


-0.018 


1 


.402 


6071 


5940 


-2.29 


3.47 






LP 815-43 


10.774 


0.272 


0.054 


0.376 


2 


.623 


0.032 


1 


.117 


6558 


6508 


-2.67 


4.17 






G 011-044 


11.091 


0.335 


0.061 


0.248 


2 


.588 


-0.009 


1 


.381 


6023 


5967 


-2.09 


4.80 






G 013-009 


9.916 


0.292 


0.054 


0.369 


2 


.609 


0.019 


1 


.195 


6402 


6319 


-2.27 


3.68 






G016-013 


9.951 


0.387 


0.129 


0.286 


2 


.583 


0.008 


1 


.627 


5676 


5529 


-0.76 


4.60 






G 018-039 


10.392 


0.346 


0.073 


0.286 


2 


.581 


-0.007 


1 


.412 


5926 


5894 


-1.52 


4.32 






G 020-008 


9.948 


0.356 


0.047 


0.251 


2 


.574 


0.001 


1 


.465 


5879 


5831 


-2.28 


4.58 


4.56 


0.43 


G 024-003 


10.393 


0.346 


0.062 


0.268 


2 


.585 


0.017 


1 


.414 


5925 


5894 


-1.67 


4.51 






G 029-023 


10.230 


0.339 


0.059 


0.332 


2 


.590 


0.008 


1 


.400 


6010 


5923 


-1.80 


3.49 






G 053-041 


11.022 


0.356 


0.083 


0.271 


2 


.589 


0.004 


1 


.460 


5848 


5810 


-1.34 


4.54 






G 064-012 


11.330 


0.282 


0.046 


0.330 


2 


.622 


0.030 


1 


.114 


6458 


6563 


-3.17 


4.54 






G 064-037 


11.144 


0.299 


0.057 


0.329 


2 


.626 


0.009 


1 


.235 


6328 


6308 


-3.12 


4.21 






G 066-030 


10.950 


0.287 


0.077 


0.354 


2 


.634 


0.018 


1 


.177 


6381 


6311 


-1.52 


4.25 






G 126-062 


9.478 


0.330 


0.063 


0.327 


2 


.588 


-0.004 


1 


.455 


6060 


5825 


-1.64 


3.75 


4.11 


0.37 


G 186-026 


10.829 


0.306 


0.041 


0.339 


2 


.608 


0.011 


1 


.229 


6280 


6266 


-2.62 


3.73 


5.20 


0.50 



difference but also with the largest uncertainty of Afy,par, 
the mean difference (A — My^-phot — -My.par) is < A >= 
—0.07 mag with a standard deviation of (t(A) ~ 0.32 mag. 
The corresponding error in My is ±0.22 mag. This induces 
an error of about ±0.1 dex in log g. Taking into account 
possible errors in the bolometric correction (adopted from 
Alonso et al. H1995|l ') and the mass as derived by interpo- 
lating in the My- log T^s diagram between the a-element 
enhanced evolutionary tracks of VandenBerg et al. H2()0()(l . 
we estimate that the error of logg is about ±0.15 dex. 

4. Abundances 

The determination of abundances is based on a-element 
enhanced ([a/Fe]= +0.4, a = O, Ne, Mg, Si, S, Ca, and Ti) 
ID model atmospheres with Toff, \ogg and [Fe/H] values 
as given in Table El and a microturbulence of 1.5 kms^^. 



The models were computed with the MARCS code using 
updated continuous opacities (Asplund et al. ll997ll and in- 
cluding UV line blanketing by millions of absorption lines. 
LTE is assumed both in constructing the models and in 
deriving abundances. 



The Uppsala abundance analysis program, 
EQWIDTH, was used to calculate theoretical equiv- 
alent widths from the models. An elemental abundance 
is determined by requiring that the calculated equivalent 
width should match the observed one. When more than 
one line of the same species were measured for a star, 
the mean value is adopted by giving equal weight to each 
line. 
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4.1. Iron 

By inspecting the blue spectra, 19 apparently unblended 
Fell lines were selected for determining the iron abun- 
dance. Differential log 5/ values of these lines were derived 
by an inverted abundance analysis of four 'standard' stars, 
HD 103723, HD 160617, G 13-09 and HD 140283 adopt- 
ing iron abundances of [Fe/H]= -0.79, -1.79, -2.30 and 
—2.42, respectively, from Nissen et al. H2002|l . In this cal- 
culation, only lines with equivalent widths between 5 and 
50 mA were included. If a differential log gf value is avail- 
able for more than one of the four standard stars, then the 
mean value (given in Table is adopted. We note that 
the [Fe/H] values of Nissen et al. 1)200211 are based on very 
weak Fe 11 lines in the red spectral region measured from 
extremely high S/N spectra and analyzed differentially to 
the solar flux spectrum. The same procedure could not be 
applied to the blue Fe 11 lines because they are quite strong 
in the solar spectrum and often blended with weak lines 
making the solar equivalent width measurement difficult. 

In calculating abundances from the Fell lines we 
adopted the Unsold (1955) approximation to the Van der 
Waals interaction constant with an enhancement factor 

— 2.5. But the adopted value of does not have a 
large effect on [Fe/H] since most Fe 11 lines in the stars are 
weak, especially in metal-poor stars. For the most metal- 
rich stars in our sample the change of [Fe/H] is around 
-fO.OSdex when E^ is changed from 2.5 to 1.5. 

In the more metal-poor stars ([Fc/H] < —1.5) the Fell 
lines are so weak that the derived mctallicity is practi- 
cally independent of the microturbulcnce. For such stars 
we have assumed ^micro — 1-5 kms^^. For the more metal- 
rich stars ^micro has been determined by requesting that 
the derived [Fe/H] values should be independent of equiv- 
alent width. The error of ^micro is about ±0.2 kms^^. 

The present sample has 10 stars in common with the 
sample of Nissen et al. ( 2002 ) . The metallicities of these 
stars range from [Fe/H] = —0.75 to —2.70. The mean 
difference (present - 2002) is < A[Fc/H] >= -0.01 with a 
standard deviation of ±0.04 dex. This small scatter shows 
that precise [Fc/H] values have been determined in both 
papers. 

4.2. Sulphur 

The log gf values of the S i lines listed in Table El are 
adopted from the Coulomb approximation calculations of 
Lambert & Luck H1978|l . From the equivalent widths of 
the AA8694.0, 8694.6 Si lines in the solar flux spectrum 
(Kurucz et. al. I1984tt of 11.0 and 28.8 mA, respectively, 
we derive a solar sulphur abundance of loge(S) — 7.20. 
This is also the meteoritic abundance (Grevesse & Sauval 
I1998|l . but given possible errors in the log gf values (the 
NIST database lists values that are 0.04 dex higher than 
those of Lambert & Luck) the exact agreement is quite 
fortuitous. Anyhow, when deriving the differential sulphur 
abundance with respect to the Sun, i.e. [S/H], a possible 
error in log gf is of no concern. Furthermore, since the 



Table 2. The list of lines used to determine the abun- 
dances of Fe, S and Zn. Measured equivalent widths 
are given for three representative stars: (a) HD 194598, 
[Fe/H] = -1.17, (b) HD 160617, [Fe/H] = -1.79, and (c) 
HD 140283, [Fe/H] = -2.42. 



ID Wavelength Exc. Pot. log gf W{a.) Wih) W{c) 





A 


eV 




mA 


mA 


mA 


Fen 


4128.74 


2.58 


-3.73 


15.0 


5.0 




Fen 


4178.86 


2.58 


-2.61 


55.8 


39.7 


17.8 


Fen 


4233.16 


2.58 


-2.01 


81.3 


68.1 


41.8 


Fen 


4416.83 


2.78 


-2.65 


46.3 


29.5 


11.2 


Fen 


4489.18 


2.83 


-2.96 


32.5 


17.3 


5.5 


Fen 


4491.41 


2.85 


-2.80 


37.9 


21.4 


7.4 


Fen 


4508.29 


2.85 


-2.41 


55.1 


37.4 


16.2 


Fen 


4515.34 


2.84 


-2.56 


48.9 


31.2 


12.4 


Fen 


4520.23 


2.81 


-2.66 


46.0 


28.9 


11.0 


Fen 


4522.63 


2.84 


-2.22 


64.6 


48.3 


22.5 


Fen 


4541.52 


2.85 


-3.04 


28.4 


14.4 


4.5 


Fen 


4555.89 


2.83 


-2.43 


53.9 


37.5 


16.1 


Fen 


4576.34 


2.84 


-3.01 


28.5 


15.0 


5.0 


Fen 


4582.83 


2.84 


-3.24 


20.9 


9.8 


3.2 


Fen 


4583.83 


2.81 


-1.91 


82.1 


63.8 


37.9 


Fen 


4620.52 


2.83 


-3.36 


17.5 


7.4 


2.5 


Fen 


4656.98 


2.89 


-3.66 


8.1 


3.1 




Fen 


4666.75 


2.83 


-3.31 


17.9 


7.5 


2.9 


Fen 


4923.93 


2.89 


-1.45 




82.9 


58.5 


Znl 


4722.16 


4.03 


-0.39 


23.8 


8.6 


3.9 


Znl 


4810.54 


4.08 


-0.17 


29.0 


13.5 


5.5 


Si 


8693.96 


7.87 


-0.56 


1.9 






Si 


8694.64 


7.87 


0.03 


5.1 


4.1 




Si 


9212.87 


6.52 


0.38 


77.0 


57.4 


19.8 


Si 


9237.54 


6.52 


0.01 


53.2 


34.4 


9.2 



AA8694.0, 8694.6 Si lines are quite weak and apparently 
without any blends in the solar spectrum, the derived 
[S/H] is not subject to the kind of systematic errors which 
may affect [Zn/H] as discussed in Sect 14.31 

The flux around the A9237.5 Si line is slightly de- 
pressed relative to the true continuum by the broad wing 
of the Paschen C Hi line at 9229 A. In this connection, we 
note that the equivalent widths of the S i lines were mea- 
sured relative to the local, apparent continuum around 
the lines. Furthermore, when deriving the sulphur abun- 
dance we neglected the contribution of the H i line to the 
line absorption coefficient. For a few representative stars 
a spectrum synthesis of the region around the A9237.5 S i 
line was carried out including the line absorption contribu- 
tion from the Paschen C H i line calculated as described by 
Seaton H1990|l . This exercise shows that the correction to 
the sulphur abundance derived from the equivalent width 
of the A9237.5 Si line is at most about +0.07 dex for our 
hottest stars and decreases with decreasing Tcff. In the 
case of the A9212.9 S i line, which is further away from the 
Hi line, the corresponding correction is negligible. An em- 
pirical confirmation of the effect is seen in Fig.|Sl where the 
difference of sulphur abundances derived from the A9212.9 
and A9237.5 fines is plotted as a function of Teg. Given 
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Table 3. The derived values of the effective temperature, 
surface gravity, microturbulence, and the abundances of 
iron, sulphur and zinc. 



5500 



5800 



5100 
n„(K) 



6400 



6700 



Fig. 5. The difference in S abundances derived from the 
A9212.9 and A9237.5 Si hues. The least squares fit to the 
data is shown. The rms deviation between the two sets of 
abundances is ±0.06 dex. 

the small size of the effect and the uncertainty in the line 
broadening theory of this high-series Paschcn Hi line, we 
have not included any correction for the H i line absorp- 
tion. In this connection, we note that for the two most 
metal-poor stars in our sample, G 64-12 and G 64-37, only 
the A9212.9 S i line could be detected and used for deriving 
the S abundance. 

The AA9212.9, 9237.5 Si lines are too strong in the 
solar spectrum and the spectral region is too crowded with 
telluric H2O lines to allow a reliable solar S abundance to 
be determined from these lines. In order to check that 
this pair provides S abundances on the same scale as the 
AA8694.0, 8694.6 lines, we make use of the fact that for 18 
of the 34 program stars the S abundance has been derived 
from both the AA8694.0, 8694.6 and the AA9212. 9, 9237.5 
pairs. The mean abundance difference is 0.03 dex and the 
standard deviation is ±0.08 dex. Hence, we are confident 
that the two sets of S i lines provide [S /H] values on the 
same scale. 

4.3. Zinc 

The zinc abundance is based on two Zni lines at 4722.2 
and 4810.5 A. The loggf values (-0.39 and -0.17 respec- 
tively) were adopted from the calculations of Bicmont & 
Godefroid H1980I) . From equivalent widths measured in 
the solar flux spectrum (Kurucz et. al. .1984|J of 70.3 and 
78.3 mA, respectively, and using a solar model atmosphere 
computed with the same code as that for the stars, we de- 
rive a solar zinc abundance of loge(Zn) = 4.57. This is 0.10 
dex lower than the meteoritic abundance of 4.67 ±0.04 
(Grevesse & Sauval 19981 • The same problem was encoun- 
tered by Biemont & Godefroid t j,980J , who derived a so- 
lar photospheric zinc abundance of loge(Zn) = 4.60 from 
six Zni lines (including A4722.2 and A4810.5) using the 
Holweger-Miiller (1974) model of the Sun. This difference 
between photospheric and meteoritic Zn abundances may 
well be due to a systematic error in the log gf values of the 
Zn I lines. However, by determining differential Zn abun- 
dances with respect to the Sun, i.e. [Zn/H], using our solar 
photospheric abundance of log e(Zn) = 4.57 as a reference, 
errors in log gf will cancel. On the other hand, it should be 



Star 


Toff 


logs 


^turb 


[Fe/H] 


[S/H] 


[Zn/H] 




[K] 


[cgs] 


[km/s] 








BD- 13° 3442 


6500 


4.16 


1.5 


-2.61 


-2.31 




CD-30° 18140 


6272 


4.13 


1.5 


-1.88 


-1.63 


-1.82 


CD-35° 14849 


6125 


4.11 


1.5 


-2.41 


-2.10 


-2.27 


CD-42° 14278 


5812 


4.25 


1.5 


-2.12 


-1.82 


-2.04 


HD 103723 


6040 


4.26 


1.3 


-0.82 


-0.84 


-0.95 


HD 105004 


5919 


4.36 


1.2 


-0.86 


-0.85 


-0.86 


HD 106038 


5919 


4.30 


1.2 


-1.42 


-1.00 


-1.25 


HD 108177 


6034 


4.25 


1.5 


-1.74 


-1.39 


-1.65 


HD 110621 


5989 


3.99 


1.5 


-1.66 


-1.34 


-1.63 


HD 121004 


5595 


4.31 


1.0 


-0.77 


-0.43 


-0.58 


HD 140283 


5690 


3.69 


1.5 


-2.42 


-2.11 


-2.35 


HD 146296 


5671 


4.17 


1.2 


-0.74 


-0.54 


-0.69 


HD 148816 


5823 


4.14 


1.2 


-0.73 


-0.49 


-0.61 


HD 160617 


5931 


3.77 


1.5 


-1.79 


-1.40 


-1.81 


HD 179626 


5699 


3.92 


1.2 


-1.14 


-0.85 


-1.14 


HD 181743 


5863 


4.32 


1.5 


-1.93 


-1.59 


-1.94 


HD 188031 


6054 


4.03 


1.5 


-1.79 


-1.47 


-1.76 


HD 193901 


5672 


4.38 


1.0 


-1.12 


-0.95 


-1.18 


HD 194598 


5906 


4.25 


1.3 


-1.17 


-0.99 


-1.22 


HD 215801 


6005 


3.81 


1.5 


-2.29 


-2.01 


-2.21 


LP 815-43 


6533 


4.25 


1.5 


-2.67 


-2.43 




G 011-044 


5995 


4.29 


1.5 


-2.09 


-1.71 


-1.90 


G 013-009 


6360 


4.01 


1.5 


-2.27 


-1.81 


-2.20 


G016-013 


5602 


4.17 


1.0 


-0.76 


-0.47 


-0.81 


G 018-039 


5910 


4.09 


1.5 


-1.52 


-1.15 


-1.44 


G 020— 008 


5855 


4.16 


1.5 


—2.28 


— 1.85 


—2.15 


G 024-003 


5910 


4.16 


1.5 


-1.67 


-1.29 


-1.61 


G 029-023 


5966 


3.82 


1.5 


-1.80 


-1.50 


-1.83 


G 053-041 


5829 


4.15 


1.3 


-1.34 


-1.05 


-1.29 


G 064-012 


6511 


4.39 


1.5 


-3.17 


-2.81 




G 064-037 


6318 


4.16 


1.5 


-3.12 


-2.78 




G 066-030 


6346 


4.24 


1.5 


-1.52 


-1.28 


-1.52 


G 126-062 


5943 


3.97 


1.5 


-1.64 


-1.31 


-1.71 


G 186-026 


6273 


4.25 


1.5 


-2.62 


-2.43 





noted that the Zni lines in the solar spectrum are quite 
strong and hence subject to errors in the adopted solar 
flux microturbulence (1.15kms~^) and the Van der Waals 
damping constant (an enhancement factor E-^ — 1.5 rela- 
tive to the Unsold approximation was adopted) . The solar 
equivalent widths of the Zni AA4722. 2, 4810.5 lines are 
also quite uncertain due to blends in their wings. Hence, 
the scale of the [Zn/H] values for our metal-poor stars 
could be in error by up to 0.10 dex. 



4.4. Results 

The abundances [Fe/H], [S/H] and [Zn/H] derived for the 
stars in our sample are listed in Table |31 together with the 
stellar atmospheric parameters. For the most metal-poor 
stars the Zni lines are too weak (W < 2mA) to provide 
reliable Zn abundances, and the value of [Zn/H] is not 
given. 
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In Fig. El [S/Fe] is plotted as a function of [Fe/H]. All 
our program stars have space velocities that are typical of 
halo stars, and are shown with filled circles. In addition, 
we have plotted disk stars from Chen et al. I)2U(J2|I . who 
determined sulphur abundances from five weak S I lines 
including the AA8694.0, 8694.6 pair. 

Fig. 13 shows [Zn/Fe] vs. [Fe/H]. Here the data for the 
disk stars are from Chen et al. 1)2003)1 . who observed the 
weak Znl line (xoxc = 5.79 eV) at 6362.35 A and deter- 
mined [Zn/H] by a differential model atmosphere analysis 
with respect to the Sun. This Zn line lays in the midst of 
a very broad and weak Cai auto-ionization line (Mitchell 
& Mohler'1965(l. but its equivalent width can still be mea- 
sured reliably with respect to the local continuum. In the 

[Mil 



solar fiux spectrum (Kurucz et al. I1984|l its equivalent 
width is close to 22.0 mA and the line appears unblended. 
Hence this line is highly suitable for determining Zn abun- 
dances in disk stars. The A6362.35 line is not covered by 
the UVES spectra of the present program, but in the very 
high S/N spectra of Nissen et al. (|2002(l we were able 
to detect the line for three of the more metal-rich halo 
stars and to derive [Zn/H] in a differential analysis with 
respect to the Sun: HD 103723, W = 4.9mA, [Zn/H] = 
-0.85; HD 106038, W = 2.6mA, [Zn/H] = -1.18; and 
HD 121004, W = 9.5mA, [Zn/H] = -0.53. These Zn 
abundances compare reasonably well with those given in 
Table 01 (the differences being < 0.1 dex), indicating that 
the values of [Zn/H] derived in the halo stars from the 
Zni AA4722.2, 4810.5 lines are approximately on the same 
scale as [Zn/H] for the disk stars. 

Finally, Fig.|Hlshows [S/Zn] vs. [Zn/H]. Here the disk 
stars are those from Chen et al. (2002, 2003,1 for which 
both S and Zn abundances are available. 



4.5. Statistical abundance errors 

The main contributions to the random errors in the abun- 
dances arise from the errors in the equivalent width mea- 
surements and the uncertainties in the model atmosphere 
parameters. The level of the latter contribution has been 
estimated by changing Toff of the stellar models by 100 K, 
log 5 by 0.2 dex and the metallicity by 0.2 dex. Table 01 
shows the corresponding changes in the various abundance 
ratios. Evidently, ATcft — 100 K has a negligible effect of 
[Fe/H], and rather small effects on the other abundance 
ratios, except [S/Zn]. [Fe/H] and [S/H] are quite sensitive 
to logg but the \ogg effects nearly cancel in the S/Fe ratio. 
We also see that the effect of changing the metal content of 
the model atmosphere is negligible. Finally, we note that 
the possible error in the microturbulence parameter only 
has a significant effect on the derived abundances for the 
more metal-rich stars. 

Errors in the measured equivalent widths were dis- 
cussed in Sect. 12.31 and found to be about ±0.6 mA for 
the Fell and the Zni fines as weU as the AA8694.0, 8694.6 
S I pair, whereas the error for the AA9212.9,9237.5A Si 
lines is ±2. 3mA. Adopting these errors, and taking into 
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Fig. 6. [S/Fe] vs. [Fe/H]. Filled circles refer to halo stars 
from the present program. Open circles are disk stars with 
abundances from Chen et al. H2002|l . The error bars indi- 
cate l-CT statistical errors of [S/Fe] and [Fe/H]. 
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Fig. 7. [Zn/Fe] vs. [Fe/H]. Filled circles: halo stars; open 
circles: disk stars from Chen et al. (2003). The error bars 
indicate 1-cr statistical errors of [Zn/Fe] and [Fe/H]. 



account the number of lines observed for a given star and 
the contribution from the estimated errors in T^s (±70 K), 
\ogg (±0.15 dex) and ^turb (±0.2kms^^), we have esti- 
mated individual abundance errors for each star. These 
errors are indicated with 1-sigma error bars in Figs. El 
13 andlHl At the lowest metallicities the dominating error 
contribution comes from the measurement of the equiva- 
lent widths of the weak absorption lines, whereas errors in 
the atmospheric parameters give the largest contribution 
for the more metal-rich stars. 



4.6. Systematic abundance errors 

In this section we discuss systematic abundance errors 
arising from the modelling of the line formation processes 
using plane parallel, homogeneous model atmospheres and 
the assumption of LTE. 
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from the identical calculations presented by Asplund & 
Garcia Perez H2001|l and correspond to the average of six 
Fell lines. The spectrum syntheses assume [S/Fe] = -1-0.4 
and [Zn/Fe] = 0.0 for [Fe/H] < -1, but the results are 
almost identical for reasonable [S/Fe] and [Zn/Fe]. It can 
be seen that the maximum 3D effect occurs for the most 
metal-poor stars: the 3D abundance is higher than the 
ID value by 0.07 to 0.13 dex in the case of the AA9212.9, 
9237.5 Si lines. Interestingly, however, there is a similar 
3D effect on the Fe abundance derived from Fell lines. 
Hence, it appears that we are in the favourable situation 
that the S/Fe ratio is quite immune to ID - 3D effects, 
with net differences of A[S/Fe] ^ 0.05 dex for the stellar 
parameters typical of our sample. Similarly, the 3D effects 
on the Zn/Fe ratio are very small (A[Zn/Fe] ^ 0.05 dex). 
The small 3D effects and the fact that they in general have 
the same sign are a consequence of the similarity in line 
formation depths of these Si, Zni and Fell lines. 



Table 4. Changes in derived abundances resulting from 
the listed changes in model atmosphere parameters. 







A[t] 


A[|] 


A[t] 


A[i] 


A[t] 


A[|r] 


ATeff = 


100 K 


0.01 


-0.04 


0.04 


-0.05 


0.03 


-0.08 


Alog g = 


0.2 


0.07 


0.06 


0.03 


-0.01 


-0.04 


0.03 


A [Fe/H] 


= 0.2 


0.01 


0.01 


0.01 


0.00 


0.00 


0.00 


A^turb = 


= 0.3km/s'' 


-0.06 


-0.02 


-0.02 


0.04 


0.04 


0.00 



" The abundance changes caused by changing ^turb by 
0.3kms~^ are maximum values valid for the more metal-rich 
stars. 



4.6.1. The influence of stellar granulation 

Much attention has recently been paid to the effect on 
abundance determinations of three-dimensional (3D) hy- 
drodynamical modelling of convection and granulation in 
stellar atmospheres. As first shown by Asplund et al. 
iPMUjl . there is a very significant surface cooling in metal- 
poor stars caused by the near adiabatic cooling of rising 
and expanding elements coupled with the lack of radiative 
heating due to the low line opacity. The effect on lines 
formed in the upper atmosphere can be very large. As an 
example, Asplund & Garcia Perez H2001(l estimated 3D ef- 
fects on oxygen abundances derived from OH lines to be of 
the order of —0.6 dex for turnoff stars with [Fe/H] < —2.0. 

To investigate the 3D effects on the sulphur and zinc 
lines used in the present study a differential study simi- 
lar to that described by Asplund & Garcia Perez (2001) 
has been carried out. The results are given in Table 
for two sets of model atmospheres with metallicities of 
[Fe/H] = 0.0,-1.0,-2.0 and -3.0. The first set has Tcff 
and log 5 close to the solar values, whereas the second 
set represents stars close to the turnoff of the halo pop- 
ulation. For each 3D model atmosphere, the S and Zn 
abundances are those reproducing the equivalent widths 
computed with ID MARCS with identical stellar param- 
eters and spectral line data; for Fe, the results are culled 



4.6.2. Non-LTE effects 

The S and Fe abundances have been derived from lines 
belonging to the main ionization stages of the elements, 
i.e. the neutral stage for S and the first ionized stage for 
Fe. Furthermore, the lines correspond to high excitation 
levels and are weak in our stars. Hence, they are formed 
deep in the stellar atmospheres where only small depar- 
tures from LTE are expected. This has been confirmed in 
the case of the Fell lines by Thevenin & Idiart (|1999fl . 
According to Takada-Hidai et al. (|2002|l non-LTE effects 
on the weak, high excitation Si lines at 8694 A are es- 
timated to be small (< 0.05 dex). The good agreement 
found in Sect 14.21 between sulphur abundances derived 
from the AA8694.0, 8694.6 pair and the AA9212.9, 9237.5 
lines suggests that non-LTE effects on the 9212-9238 Si 
triplet are also small. 

The case of zinc is more difficult. With an ionization 
potential, Xion(Zni) = 9.39 eV, there are approximately 
equal numbers of neutral and ionized zinc atoms at the 
temperatures and electron pressures of the line forming 
regions in the atmospheres of our stars. Hence, an over- 
ionization of Zni relative to LTE, as is the case for Fei 
(Thevenin & Idiart 19991, would lead to an underestimate 
of the abundance of Zn in the LTE analysis of Zni lines. 
However, such an over-ionization effect on Zni is likely 
to be smaller than the equivalent effect for Fei, since 
A^Zni ^ Nznii whereas A^fci << A^Fch- Furthermore, there 
may be departures from LTE in the excitation balance of 
Zni, since the lower level from which the lines originate 
is the first excited level of Zni. The discrepancy noted 
above between the meteoritic Zn abundance loge(Zn) — 
4.67 and the photospheric value of 4.57 may be due to such 
non-LTE problems. Given the complex structure of the Zn 
atom, a thorough study of non-LTE effects on the deter- 
mination of Zn abundances is clearly needed but, to our 
knowledge, none has been carried out yet. Here we simply 
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Table 5. The effects of 3D hydrodynamical model atmospheres on the derived sulphur, zinc and iron abundances 
relative to those estimated with classical ID models (A(loge) — loge(3D) — loge(lD)). Both the ID and 3D abundances 
have been computed assuming LTE. It should be noted that the Si 9212.9 A and Si 9237.5 A lines are very strong 
{W > 100mA ) at [Fe/H] = 0.0 and hence are sensitive to the spectral line broadening. Similarly, the Si 8694.6 A line 
is very weak (W < 0.5 mA ) at [Fe/H] = -3.0. 



r,ff logs [Fe/H] 
[K] 



A(loge) = loge(3D) - log£(lD) 



Si 8694.6 A Si 9212.9 A Si 9237.5 A Zni 4722.2A Zn 14810.5 A Feii 



5767 


4.44 


-0.0 


5822 


4.44 


-1.0 


5837 


4.44 


-2.0 


5890 


4.44 


-3.0 


6191 


4.04 


-0.0 


6180 


4.04 


-1.0 


6178 


4.04 


-2.0 


6205 


4.04 


-3.0 



-0.04 -0.01 

+0.02 +0.02 

+0.06 +0.13 

+0.05 +0.08 

+0.03 +0.11 

+0.04 +0.01 

+0.04 +0.09 

+0.04 +0.08 



-0.01 +0.02 

+0.02 +0.01 

+0.11 +0.09 

+0.08 +0.08 

+0.09 -0.08 

+0.03 -0.02 

+0.09 +0.03 

+0.07 +0.04 



+0.02 -0.02 

+0.00 +0.04 

+0.09 +0.10 

+0.08 +0.09 

-0.06 -0.01 

-0.03 +0.02 

+0.03 +0.06 

+0.04 +0.07 



note that, by analyzing our Zni lines relative to the solar 
flux spectrum, the problem is somewhat reduced. 

Ideally, the non-LTE calculations should be done in 
combination with 3D models. This is, however, a very 
demanding task, which has only recently become fea- 
sible (e.g. Kiselman & Nordlund 119951 Kiselman 119971 
Uitenbroek [r998 Asplund et al. 12003(1 . For complex atoms 
like Zn, the computationally less demanding 1.5D non- 
LTE problem'^ would be a good starting point, as recently 
attempted for Fel (Shchukina & Bueno [?001|l . It is note- 
worthy that the 3D non-LTE results can be significantly 
different from both the ID non-LTE and the 3D LTE re- 
sults. There are, however, very good reasons why quite 
small ( ^ 0.1 dex) 3D non-LTE corrections are expected 
in the case of S i. The S i lines are formed in deep atmo- 
spheric layers where the differences between ID and 3D 
model atmospheres are small. One therefore expects sim- 
ilar non-LTE abundance corrections in ID and 3D, as for 
example appears to be the case for O i. 

5. Discussion 

5.1. Sulphur as an a-capture element 

As can be seen from Fig. |S1 the halo stars are distributed 
around [S/Fe] ~ +0.3 except for two deviating stars, 
HD 103723 and HD 105004, that are known to have solar- 
like [a/Fe] ratios as discussed in the next paragraph. Using 

^ The approximation of treating an ensemble of vertical at- 
mospheric columns extracted from a 3D model as separate 
plane-parallel ID model atmospheres on which the ID non- 
LTE line formation calculations are performed before averag- 
ing, is normally referred to as a 1.5D non-LTE calculation. Such 
a procedure ignores all horizontal radiative transfer in the 3D 
model atmosphere. While the results are not quantitatively ex- 
act, they can nevertheless provide a qualitative assessment of 
the general behaviour of the full 3D non-LTE problem and 
yield a useful estimate of the 3D non-LTE effects. 



a maximum likelihood program that takes into account in- 
dividual errors in both x and y, we obtain a fit for the halo 
stars (excluding HD 103723 and HD 105004) 

[S/Fe] = 0.240 (±0.024) - 0.042 (±0.014) • [Fc/H] 

with a reduced chi-square xfed — If HD 103723 and 
HD 105004 are included the relation becomes 

[S/Fe] 0.153 (±0.023) - 0.085 (±0.013) • [Fc/H], 

but then Xred increases to an unacceptable high value of 
2.3. In any case the slope is much smaller than the corre- 
sponding slope of the fit to the disk star data: 

[S/Fe] = -0.014 (±0.015) - 0.308 (±0.040) • [Fe/H]. 

Hence, there is no tendency in our halo star data for a con- 
tinuing strong linear rise of [S/Fe] with decreasing [Fe/H], 
as claimed by Israelian & Rebolo 1(2001(1 . At [Fe/H] = -2 
their data suggest [S/Fe] ~ +0.7, whereas our relation cor- 
responds to [S/Fe] = +0.32. We note in this connection 
that their sample includes only three stars with [Fc/H] < 
— 1.8 and none below [Fe/H] — —2.3, whereas our sam- 
ple includes 14 stars in the range —3.2 < [Fe/H] < —1.8 
and seven below [Fe/H] = —2.3. Furthermore, the er- 
ror bars on their sulphur abundance determinations are 
large due to the weakness of the Si 8694.6 A line; in our 
study we circumvent this problem with the 9212.9 and 
9237.5 A lines which are well suited for sulphur abundance 
determinations below [Fe/H] = —2. Hence, our conclu- 
sion is that there is no need to invoke element produc- 
tion by hypernovae or very massive supernovae to explain 
the general behaviour of sulphur. Our data are in good 
agreement with traditional Galactic evolution models with 
near-instantaneous production of a-elements by Type II 
supernovae and delayed production of the iron-peak ele- 
ments (e.g. Chiappini et al. 119991 Goswami & Prantzos 
EDOOj). 
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All stars from the present investigation (plotted with 
filled circles) have halo kinematics, including six stars 
with -0.9 < [Fe/H] < -0.7 which have Galactic ro- 
tational velocities of less than 50kms~^ i.e. well below 
the characteristic rotational velocity of thin disk stars 
(225kms~^) and thick disk stars (175 kms""'^). As shown 
by Nissen & Schuster H1997|l . there is an overlap be- 
tween halo and thick disk stars in the metallicity range 
— 1.0 < [Fe/H] < —0.6. Four of our six stars in this metal- 
licity range have enhanced S/Fe ratios like the thick disk 
stars, but two (HD 103723 and HD 105004) have a solar 
S/Fe ratio. These two stars also show solar a/Fe ratios in 
other a-elements, such as O, Mg, Si, Ca, and Ti (Nissen & 
Schuster I1997|l . On the basis of the stars' Galactic orbits, 
Nissen & Schuster suggested that they may have been ac- 
creted from dwarf galaxies with a chemical evolution that 
has proceeded more slowly than in the inner part of our 
Galaxy, where the 'normal' halo stars formed. The fact 
that S and the classical a-elements, Mg, Si and Ca, ex- 
hibit the same behaviour in these 'anomalous' stars is a 
further indication that sulphur belongs to this group of el- 
ements. Thus, there seems to be little ground, on the basis 
of our data, for the reservations expressed by Prochaska 
et al. C2UU0 ) concerning the use of S as an a-element in 
the analysis of abundance ratios in DLAs. 

5.2. Zinc as a tracer of iron 

Zinc is an interesting element with a number of possible 
nucleosynthesis channels: neutron capture (s-processing) 
in low and intermediate mass stars as well as explosive 
burning in Type II and la SN (Matteucci et al. HM^ . 
Furthermore, zinc is a key element in studies of elemental 
abundances of damped Lya systems, because it is one of 
the few elements which in the interstellar medium are not 
depleted onto dust and measurements of its interstellar 
absorption lines have several other practical advantages 
(Pettini, Boksenberg & Hunstead ll990|l . 

Interstellar iron, on the other hand, can exhibit large 
gas depletions, and it has thus become customary to use 
Zn as a tracer of Fe in studies aimed at investigating the 
chemical evolution and dust content of galaxies, particu- 
larly at high redshifts. The underlying assumption is that 
the abundances of Zn and Fe vary in lockstep, as orig- 
inally found by Sneden, Gratton, & Crocker H1991|) al- 
though their data, obtained with 2-3 m class telescopes, 
exhibited considerable scatter. 

Our ignorance of the nucleosynthetic origin of Zn has 
prompted some to question the validity of using it as a 
proxy for Fe. We can now reassess the validity of this as- 
sumption with modern data such as those presented here. 
Inspection of Fig.Qshows that, within the errors, [Zn/Fe] 
is indeed approximately solar, at least between [Fe/H] — 
and —2. Over this range the mean of the 61 measure- 
ments shown in Fig.dis [Zn/Fe] = +0.03 ± 0.08 (Icr). A 
similar conclusion was recently reached by Mishenina et 
al. (j2002.1 who published a survey of Zn abundances in 



90 disk and halo stars based on the equivalent widths of 
the AA 4722.2, 4810.5, 6362.35 Zni lines in high resolution 
spectra of dwarf and giant stars and concluded that their 
data "confirm the well-known fact that the ratio [Zn/Fe] 
is almost solar at all metallicities" . 

Looking more closely at the data in Fig. [T] there may 
be hints of subtle trends which, if confirmed by larger 
samples, could provide clues to the nucleosynthetic ori- 
gin of Zn. For example, Prochaska et al. I|2000|l claimed 
that in thick disk stars [Zn/Fe] ~ -1-0.1. This mild over- 
abundance of Zn relative to Fe in thick disk stars is also 
present in the Mishenina et al. H2002|l sample, once the 
stars are separated on the basis of their kinematics into 
halo, thick and thin disk populations, as recently pointed 
out by Nissen l)2()()3|l . Furthermore, there seems to be a 
gradient in [Zn/Fe] as a function of [Fe/H] for the halo 
stars in the Mishenina et al. sample with the highest values 
of [Zn/Fe] being measured in the most metal-poor stars. 
The few stars in the present study with [Fe/H] < —2.0 
show a similar effect, with a mean [Zn/Fe] ~ +0.1 (see 
Fig. [7||. On the other hand, before making too much of 
these trends, it is important to bear in mind that there 
may well be systematic errors at the ±0.10 dex level in 
the determinations of [Zn/Fe], as discussed in Sect. 14.31 
Clearly, [Zn/Fe] in halo and disk stars should be stud- 
ied further, both observationally and in terms of Galactic 
chemical evolution models. 

Notwithstanding these uncertainties, it can be seen 
from Fig.lHlthat the [S/Zn] ratio shows the same trend vs. 
metallicity as the classical [a/Fe] trend, i.e. a plateau at 
[S/Zn] ~ +0.25 until [Zn/H] ~ -1.0 and then a decline of 
[S/Zn] to the solar ratio at [Zn/H] ~ 0. Complications may 
arise, however, in the transition region between the halo 
and the disk where the 'low'-a stars identified to date oc- 
cur. More measurements in this metallicity regime would 
be highly desirable. 

5.3. Sulphur and zinc in damped Lya systems 

Armed with the results of Fig. |H1 we are now in a posi- 
tion to examine the [S/Zn] ratio in DLAs and compare 
it to the values measured in Galactic stars. Despite the 
large database of abundance measurements in DLAs ac- 
cumulated over the last few years, there are still rela- 
tively few determinations of [S/Zn] in DLAs. The reasons 
are two-fold. First, the rest wavelengths of the Sii triplet 
AA1250, 1253, 1259 are so close to the wavelength of Lya 
(1215.67 A) that the S ii lines often fall within the Lya for- 
est, where blending can be a problem. Second, with rest 
wavelengths AA2026,2062, the Znii doublet lines are sep- 
arated by ~ (1 + Zabs) X 800 A from the Sii triplet; thus, 
at the redshifts Zabs = 2 — 3 of most DLAs, the two sets of 
spectral features fall in widely separated regions of the op- 
tical spectrum which may not be covered simultaneously 
by some instruments. For this reason, most measurements 
of the [S/Zn] ratio in DLAs have become available only 
recently thanks to the wide spectral coverage and high ef- 



P.E. Nissen et al.: Sulphur and zinc abundances in Galactic stars and damped Lya systems 
Table 6. Interstellar sulphur and zinc measurements in damped Lya systems. 
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logiV(Znii) 
(cm- 2) 


[Zn/H]" 


[S/Zn]'' 


Ref.'^ 












yUUUU— 2d2U 


0.39Ui 


oi A 1 _L n no 


14. (0 ± U.Uo 


1 o m _L n nK 
12. Ui ± U.Uo 


o n'7 
— 2.0/ 


±U.lD 


1 O 

1,2 


I^UUio — UU4 


i.y ( zyb 


2U.80 ± U.Uo 


1 OQ _i_ n no 
io.zo ± U.Uo 


1 o 7/1 _i_ n (\A 
12. ( 4 zh U.U4 


— U. 10 


1 n ni 
+U.U1 


Q 
O 


/^m nn i i on 




Zi.4U ± U.Uo 


1 c 1 o _L n nc 
io.io ± U.Uo 


12. 4o ± U.IU 


— l.t)2 


±U.l0 
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Q0407-4410 


2.5505 


21.13 ±0.10 


14.82 ±0.06 


12.44 ± 0.05 


-1.36 


-0.15 


5 


Q0407-4410 


2.5950 


21.09 ±0.10 


15.19 ±0.05 


12.68 ±0.02 


-1.08 


-0.02 


5 


Q0551-366 


1.96221 


20.50 ±0.08 


15.38 ±0.11 


13.02 ±0.05 


-0.15 


-0.17 


6 


0841+1256 


2.3745 


21.00 ±0.10 


14.77 ±0.03 


12.20 ±0.05 


-1.47 


±0.04 


7 


HE2243-6031 


2.33000 


20.67 ±0.02 


14.88 ±0.01 


12.22 ±0.03 


-1.12 


±0.13 


8 


B2314-409 


1.8573 


20.90 ±0.10 


15.10 ±0.05 


12.52 ±0.03 


-1.05 


±0.05 


9 


Q2343+125 


2.4313 


20.35 ±0.05 


14.71 ±0.08 


12.45 ±0.06 


-0.57 


-0.27 


3,10 



^ [log (Zn/H)DLA - log (Zn/H)Q], where log (Zn/H)Q = -7.33 (Grevesse & Sau vallT^ . 

[log (S/Zn)DLA - log (S/Zn)0], where log (S/Zn)© = 2.53 (Grevesse & Sauval lT?^ . 

" References— 1: Molaro et al. (;20'0nt : 2: Lu, Sargent, & Barlow (1998); 3: Petitjean, Srianand & Ledoux 1120021 
4: Prochaska et al. 1200111 : 5: Lopez & Ellison 1200311 : 6: Ledoux, Srianand, & Petitjean 120021 : 7: Centurion et 
al. 120031 : 8: Lopez et al. 120021 : 9: Ellison & Lopez 120011 : 10: Dessauges-Zavadsky, Prochaska, & D'Odorico 



ficiency at blue and red wavelengths of the VLT/UVES 
combination. 

In TableElwe have collated from the literature all mea- 
surements of [S/Zn] in DLAs from spectra obtained with 
8-10 m class telescopes; references to the original works are 
given in the last column of the Table. The total sample 
consists of ten DLAs at redshifts Zabs = 1-86 — 3.39; their 
values of [S/Zn] vs. [Zn/H] are compared with the stel- 
lar measurements in Fig. |^ The behaviour of the [S/Zn] 
ratio in DLAs is evidently not the same as that seen 
in Galactic stars. Taken together, the ten DLAs consid- 
ered here do not show evidence for an a-element enhance- 
ment; the mean and standard deviation for the sample are 
[S/Zn] = -0.01 ± 0.15. Considering only the seven DLAs 
with [Zn/H] < -1, we find [S/Zn] = ±0.05 ±0.11. 

It is hard to identify systematic effects which may 
be the cause of this apparent offset between DLAs and 
Galactic halo stars. The column density errors listed in 
Table El which in turn translate to the typical abundance 
errors illustrated in Fig.^ are those quoted by the authors 
of the original papers, as referenced in the Table. These 
errors are the Icr uncertainties returned by the absorption 
line fitting computer codes and generally reflect the ran- 
dom errors in the line equivalent widths. They are prob- 
ably underestimates of the true uncertainties (sec, for ex- 
ample, the discussion of this point by Kirkman et al. 2003.) , 
but this should result in an increased scatter of the data 
points in Fig.^ rather than a systematic offset. Generally, 
the S II triplet lines are stronger than the Znii doublet, so 
that saturation may be an issue. However, with three ab- 
sorption lines available, it is normally possible to assess the 
degree of saturation reliably, unless the distribution of ab- 
sorber properties is markedly irregular (Jenkins ll986|l ; the 
line fitting programs used in the original analyses of the 
data in Table are well up to this task. Similar considera- 



tions apply to resolving the Znii AA2026.1, 2062.7 doublet 
lines from the nearby Mgi A2026.5 and Cm A2062.2 lines. 

Turning from ion column densities to element abun- 
dances, it is also difficult to find reasons why the abun- 
dance of S should have been systematically underesti- 
mated, or that of Zn overestimated. Both S ii and Zn ii 
are the major ionization stages of their respective ele- 
ments in H I regions, and corrections for unobserved ion 
stages are expected to be unimportant (e.g. Vladilo et al. 
l2UUlf) . Neither S nor Zn show much affinity for dust and 
the problem is further lessened in DLAs which generally 
show only mild depletions of even the refractory elements 
(Pettini et al. I1997|l . Errors in the solar abundance scale 
do not seem a plausible explanation for the difference ei- 
ther. Recall that while the stellar abundances are derived 
differentially relative to the Sun, those in DLAs have to 
be referred to a solar scale; here we have adopted the me- 
teoritic abundances of S and Zn from the compilation by 
Grevesse & Sauval H1998fl . However, as discussed earlier 
(Sect. 4.2 and 4.3), the Si lines used in the present study 
give a solar photospheric abundance of S which is the same 
as the meteoritic one, and while for Zn there is a O.lOdex 
offset, this offset is in the wrong direction for reconciling 
stellar and DLA data. 

In summary, on the basis of our current knowledge of 
S and Zn in the local interstellar medium, the values given 
in Table ini should reflect the true interstellar abundances 
of these two elements in the high redshift galaxies giving 
rise to the damped Lya systems. Thus, we are led to con- 
clude that the difference between Milky Way stars and 
DLAs is probably real, confirming the earlier analysis by 
Centurion et al. (2000) which was based on noisier 4-m 
data. A similar conclusion has also been reached by sev- 
eral studies which targeted depleted elements (chiefiy Si 
and Fe) and then attempted to correct for the fractions in 
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Fig. 9. [S/Zn] vs. [Zn/H] in Milky Way stars and damped 
Lya systems. Typical errors are shown. 



dust (e.g. Vladilo [20021 and references therein), although 
the method used here is clearly more direct, as it does not 
rely on the accuracy of the dust corrections. 

Interpreted within standard chemical evolution mod- 
els (e.g. Calura, Matteucci, & Vladilo .2003)) . the lack of 
a-element enhancement in DLAs may be taken as evi- 
dence for low and intermittent rates of star formation over 
their past history (relative to the time when we observe 
them). Put simply, in this scenario the overall metallic- 
ity grows only slowly with time and the iron-peak ele- 
ments released by Type la supernovae can 'catch-up' with 
the overall chemical enrichment during quiescent periods 
between isolated bursts of star formation. Such a mode 
of star formation is seen locally in dwarf irregular and 
dwarf spheroidal galaxies whose stars also do not show 
an a-element enhancement at low metallicities (Carigi, 
Hernandez, & Gilmore HUU^ Venn et al. Shetrone 
et al.EnnSI Tolstoy et al. lM7H|l . and is indeed envisaged 
in theoretical models of DLA galaxies (e.g. Mo, Mao, & 
White 1998.1 . 

Nevertheless, it is somewhat remarkable that we have 
not found even a single DLA with a high [S/Zn] ratio in 
our (albeit small) sample. At lower redshifts {z < 1), DLA 
galaxies are a heterogeneous group, exhibiting a variety of 
morphologies and surface brightnesses (Boissier, Peroux, 
& Pettini .2003) : if this is also the case at high z, we may 
have expected a wider range of [S/Zn] values. Possibly the 
issue is still clouded by small number statistics, or there 
may be more fundamental differences between high- and 
low-redshift DLAs. The findings by Kanekar & Chengalur 
H2003|l may be relevant here. These authors derived esti- 
mates of the spin temperature Tg in 24 DLAs, 11 of which 
(all at ^abs < 1) have optical identifications. Kanekar & 
Chengalur l|2()(),'-{(l find that all DLAs with high values of 
spin temperature (Tg ^ 1000 K) are identified with dwarf 
or low surface brightness galaxies, while DLAs with low 
values of Tg are invariably associated with large, luminous 
galaxies. Furthermore, low z DLAs exhibit both high and 
low values of Tg, while high redshift (zabs ^ 2) DLAs have 
preferentially high spin temperatures. The lack of a sig- 
nificant a-element enhancement in the DLAs considered 



here can be understood within the picture put forward 
by Kanekar & Chengalur ('2003 ) , if dwarf and low surface 
brightness galaxies dominate the cross-section for DLA 
absorption at high redshift, and if these types of galaxies 
generally tend to have low and intermittent rates of star 
formation. 

Centurion et al. l^IUUU\i and Vladilo (|2002|l suggested 
that there may be a mild trend of decreasing [S/Zn] 
(and more generally [a/Fe]) with increasing metallicity 
in DLAs. While the DLA data shown in Fig. are not 
inconsistent with such a possibility, the number of reliable 

5 and Zn abundance measurements needs to be consider- 
ably larger than the current sample before the reality of 
such a trend can be assessed statistically. 

6. Summary and conclusions 

One of the motivations for carrying out the present sur- 
vey of S and Zn abundances in Galactic stars is the im- 
portant role played by these two elements in deciphering 
the chemical enrichment and star formation histories of 
damped Lya systems. Up to now measurements of the 
sulphur abundance in Galactic stars have been sparse and 
doubts have been raised as to whether S is a typical a- 
capture element formed in Type II supernovae (Israelian 

6 Rebolo nnnH Takada-nidai et al. IM)^ . The lack of 
data for sulphur, coupled with insufficient knowledge of 
the nucleosynthetic channels for the production of Zn, has 
led some (e.g. Prochaska et al. l2000|l to question the valid- 
ity of using these elements, and in particular their ratio, 
in the interpretation of the abundance patterns seen in 
DLAs. 

The new data presented here have clarified the situ- 
ation as regards Milky Way stars. By targeting the S i 
AA9212.9, 9228.1, 9237.5 triplet we have overcome the lim- 
itations of most earlier studies and probed the abundance 
of sulphur with higher precision and to lower metallici- 
ties than had been possible previously. We find that the 
trend of [S/Fe] as a function of [Fe/H] is very similar to 
those of other typical a-capture elements, Mg, Si and Ca. 
[S/Fe] is nearly constant at a level of [S/Fe] ~ -|-0.3dex 
in the metallicity range —3.2 < [Fe/H] < —1.2, starts to 
decrease at [Fe/H] ~ —0.7, and reaches a solar ratio at 
[Fe/H] ^ 0. In the halo-disk transition region, at metal- 
licities — 1.2 < [Fe/H] < —0.7, there is significant scatter 
in the values of [S/Fe]. More stars in this interval should 
be studied to look for possible correlations between stel- 
lar kinematics and [S/Fe]. Precise abundance analyses by 
Fuhrmann 1,1998;) : Gratton et al. (2000J; and Fehzing et 
al. H2003|l have shown a clear separation in [Mg/Fe] be- 
tween thin and thick disk stars and it will be important 
to establish if this is also the case for [S /Fe] . 

As far as zinc is concerned, we confirm the results of 
several earlier surveys which have shown that, to a first ap- 
proximation, Zn tracks Fe over three orders of magnitude 
in [Fe/H]. There is some evidence for a small overabun- 
dance of Zn ([Zn/Fe] ~ -1-0.10) for metal-poor disk stars 
and halo stars with [Fe/H] < —2. However, there may be 
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systematic errors in [Zn/Fe] at a level of ±0.1 dex due to 
the difficulty in analyzing the AA4722. 2, 4810.5 Znl lines 
in the solar spectrum. Non-LTE effects on the derived zinc 
abundances are also a potential problem that should be 
addressed. New studies of the nucleosynthesis of zinc in 
supernovae and by the s-process in AGB stars would help 
in understanding the reasons why Zn behaves like an Fe- 
peak element and the origin of the 0.1 dex offset at low 
metallicities. 

When we turn to damped Lya systems, however, the 
Galactic pattern of S and Zn abundances does not seem to 
apply. The sample is still small (accurate measurements of 
[S/Zn] are available for only ten DLAs) and there is scat- 
ter in the data but, taken at face value, there is no obvious 
a-element enhancement in DLAs. We can not identify any 
systematic effect in the data nor in their analysis which 
would mask such an enhancement, if it were there, and 
conclude that its absence is probably real. Presumably, 
it is indicative of a bursting history of star formation in 
the galaxies giving rise to damped systems. However, it is 
important to realise that a variety of recent observations 
suggest that the a-enhancement exhibited by metal-poor 
stars of the Milky Way may in fact be the exception rather 
than the rule — it is generally not seen in dwarf spheroidal 
and dwarf irregular galaxies (Venn et al. I20U3I Shetrone 
et al.Enni Tolstoy et al. l^nUSI Aloisi et al HUU^ . in old 
stars of the Large Magellanic Cloud (HiU et al. ^10(7), 
in DLAs (this paper), in some Galactic halo stars with 
large orbits (Nissen & Schuster ll997|l . and in the globular 
cluster Pal 12 which may originally have been part of the 
SgrdSph galaxy (Cohen (2003,). The challenge now is to 
incorporate this rapidly growing set of abundance mea- 
surements into a comprehensive picture of the chemical 
evolution of galaxies. 
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